Determination of ^B(p,7)^C reaction rate from breakup 
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The astrophysical factor of the *B(p,7)^C at zero energy, 518(0), is determined from three-body model anal- 
ysis of breakup processes. The elastic breakup ■^''*Pb(^C,p*B)^''*Pb at 65 MeV/nucleon and the one-proton 
removal reaction of at 285 MeV/nucleon on C and Al targets are calculated with the continuum-discretized 
coupled-channels method (CDCC) and the eikonal reaction theory (ERT), respectively. The asymptotic normal- 
ization coefficient (ANC) of '"'C in the p-*B configuration, CpSg, extracted from the two reactions show good 
consistency, in contrast to in the previous studies. As a result of the present analysis, 5*18(0) = 66 ± 10 eVb is 
obtained. 

PACS numbers: 24.10.Eq, 25.60.-t, 21.10.Jx, 26.20.Cd 



Introduction. In low-metallicity supermassive stars, the 
proton capture reaction of ^B, ®B(p,7)^C ignites the explosive 
hydrogen burning ^ i 

«B(p, 7)9C(«,p)12n(p, 7)"0(/3+i.)"N(p, 7)^0. 

This process called hot pp chain is expected to be a possible al- 
ternative path to the synthesis of the CNO elements. Because 
of the difficulties in measuring the '*B(p,7)^C cross section 
ap-y at very low energies, several alternative reactions have 
been proposed |2-4] to indirectly determine the astrophysical 
factor (e) 

S'i8(£) = a-p^£exp[27r?7]. (1) 

Here, e is the relative energy between p and ^B in the center- 
of-mass (cm.) frame and 77 is the Sommerfeld parameter Be- 
cause an astrophysical factor has quite weak energy depen- 
dence, several previous studies have paid special attention to 
the evaluation of 5*18 (e) at zero energy, 5i8(0) 1 1- 5]. 

The Coulomb dissociation method |4] is based on the as- 
sumption that elastic breakup of by a heavy target, e.g., 
^•^^^Pb, is essentially a one-step electric dipole (El) transition 
to the p + ^B continuum. Then ap^ can be obtained by eval- 
uating the cross section of the inverse process of the breakup 
reaction |6]. This assumption needs to be examined, since nu- 
clear breakup. Coulomb dissociation with higher multipolar- 
ities, and multi-step transitions can play non-negligible roles 
even in El -dominated breakup processes fT]. In fact, an at- 
tempt to evaluate these higher-order contributions was made 
in Ref. we will return to this point later 

The asymptotic normalization coefficient (ANC) 
method f^, which is one of the most important tech- 
niques of indirect measurements have been used in several 
studies l?!, ,9.- .- 15.1 in order to determine astrophysical reaction 
rates. The basic idea of the ANC method is that only 
the tail of the overlap between the initial and final states 
contributes to a reaction at stellar energies. Thus, the purpose 
in the present case is to determine the ANC C of the ^C 
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wave function in the p + ^B configuration by using some 
alternative reactions. In Refs. |0] and llH, respectively, the 
d(^B, ^C)n reaction at 1 1.4 Mev/nucleon and the one-proton 
removal reaction of ^C at 285 MeV/nucleon were analyzed 
to determine C^fg, and hence 518(0). One of the most 
important conditions for the ANC method is that a reaction 
used to determine the ANC must be peiipheral. From this 
aspect, transfer reactions at low incident energies ll9l-ll4ll and 
nucleon removal reactions in wide range of energies f 151 have 
been used as alternative reactions for the ANC method. In 
Ref. lfl6ll . it was demonstrated that an ANC can be extracted 
from an elastic breakup cross section (angular distribution) 
for which El breakup plays a dominant role. Later this 
method was carefully examined and justified [7]; important 
findings of the work are i) El -dominated breakup processes 
are peiipheral with respect to the relative coordinate between 
the two fragments after the breakup, ii) the breakup cross 
section in a coupled-channel framework is proportional to 
the square of the ANC to be determined, and iii) if the two 
fragments are ejected in forward angles, which is the case 
in usual breakup experiments of unstable nuclei, dynamical 
excitation of each fragment during the breakup process has 
no essential effects on the ANC. 

TABLE I. Astrophysical factors of *B(p,7)'"*C in previous studies. 





5i8 [eVb] 


method 


Beaumel et al. [2J 


45 ± 13 


ANC (transfer) 


Trache et al. [3] 


46 ±6 


ANC (proton removal) 


Motobayashi [4] 


77± 15 


Coulomb dissociation 


Wiescher et al. [1] 


210 


shell model 


Descouvemont [5] 


72, 80 


cluster model 



We show in Table I] the 5i8(0) reported in the aforemen- 
tioned indirect measurements [2-4], together with theoretical 
evaluations ILJJ- One sees that the two theoretical values 
have a large difference of about factor of 3. Experimental re- 
sults seem to support the 5i8(0) obtained by a cluster model 
calculation |5]. There is, however, still a significant discrep- 
ancy of about 50% between the 5i8(0) obtained by Coulomb 
dissociation method |i4J and the ANC method ||2,l3|l- 
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FIG. 1 . Illustration of the three-body system. 



In this Rapid Communication, we reinvestigate the 
Coulomb dissociation |01 (elastic breakup) and the proton re- 
moval process lH of 9C by means of coupled-channel cal- 
culation with a three-body (p + *B + target) model. We 
adopt the continuum-discretized coupled-channels method 
(CDCC) lfT7l-[l9ll f or the former and the eikonal reaction the- 
ory (ERT) MmTl for the latter; we use the ANC method for 
both reactions. The main purpose of the present study is to 
show the consistency between the two values of S'i8(0) ex- 
tracted from these two types of breakup, and thereby deter- 
mine iS'i8(0) with high reliability. 

Theoretical framework. In Fig. [T] we show schematic il- 
lustration of the three-body (p + + target) system. The 
scattering between ^C and a target nucleus A is described by 
the Schrodinger equation 



- — V], + h + U{r,,rn)-E 



*(r,i?)^0, (2) 



where 'I'(r, i?) is the tree-body wave function and r (R) is 
the coordinate of (^C) relative to p (A). The reduced mass 
between ^C and A is denoted by /i and E is the total energy of 
the three-body system in the cm. frame. The internal Hamil- 
tonian of ^C is shown by h. The interaction [/ (rp , re ) is given 
by 

(3) 

where vj^^ and V^^^ are the nuclear and Coulomb interac- 
tions, respectively, between X and A; X represents a fragment 
particle of the projectile, i.e., p or ^B. Similarly, rx denotes 
the relative distance between X and A. 

In the present analysis of the elastic breakup of ^C, we solve 
Eq. (H with eikonal-CDCC (E-CDCC) |7, 22]. E-CDCC as- 
sumes eikonal approximation to the scattering wave between 
^C and A. As a result, the total wave function ^'(r, Ji) is ex- 
pressed by 



where ^dr) is the internal wave function of ^C with c the 
channel indices {i, £, S, I, m}; i > (i — 0) stands for 
the ith discretized-continuum (ground) state, and £, S, and / 
are, respectively, the orbital angular momentum, the channel 
spin, and the total angular momentum of the p and ^B system. 
m is the projection of / on the z-axis taken to be parallel to 



the incident beam; mp is the value of m in the i ncident c han- 
nel. b is the impact parameter defined by 6 = a/x^ + j/2 ^jj-jj 
R = {x, y, z) in the Cartesian representation. The use of the 
Coulomb incident wave 0^ [h, z) instead of the plane wave 
exp(isrc • R) in the eikonal approximation is one of the most 
important features of E-CDCC; Kc is the asymptotic wave- 
number vector of ^C in channel c from A. In the actual calcu- 
lation, we use an approximate asymptotic form of (f)^ {b.z). 
E-CDCC is shown to work very well for describing both the 
nuclear and Coulomb breakup processes with high accuracy 
and computational speed |7, 22]. 

The one-proton removal reaction, its stripping component 
in fact (see below), is analyzed by means of the eikonal reac- 
tion theory (ERT) ll20l I21I1 . which can calculate an inclusive 
cross section, such as a nucleon removal cross section, in the 
CDCC framework. ERT uses a formal solution (the scattering 
matrix S) to the coupled-channel equations of E-CDCC, and 
makes adiabatic approximation to only the nuclear part of S. 
Then one can obtain the most important result of ERT, i.e., the 
product form of S i20ll 



(5) 



where and show the contributions from the constituents 
b and c of the projectile, respectively. At this stage, how- 
ever, this result can be derived only when b or c is charge- 
less, which is not the case for the ^C projectile consisting 
of p and ^B. Therefore, in the present study, we neglect the 

Coulomb breakup process in the one-proton removal process 

(c) 

and replace the Coulomb interaction Vp (rp) with 



(6) 



Then we can calculate the one-proton removal cross section 



(T_p with 



CTbu 



(7) 



as in Refs. 120112 ih . In Eq. (|7]l, (Tbu and asti denote the elas- 
tic breakup cross section and the stripping cross section, re- 
spectively; ERT is used to evaluate asu- The accuracy of the 
replacement of Eq. (|6|l can be examined by calculating u-p 
with and without the Coulomb breakup. It is confirmed that 
the Coulomb breakup contributes to u-p for C and Al targets 
by about 5%. Thus, we conclude that the Coulomb breakup 
by these two targets can be neglected with 5% errors. Below 
we include this amount in uncertainties of 5*18 (0) extracted 
from (T_p. 

Model setting. For both the elastic breakup and one-proton 
removal processes, the p-^B wave function is calculated with 
the same Hamiltonian h. We include only the intrinsic spin of 
p. We adopt the standard Woods-Saxon central potential with 
the radial parameter Rq — 1.25 x 8^/"^ fm and the diffuseness 
parameter oq = 0.65 fm. The Coulomb interaction between 
a point charge (p) and a uniformly charged sphere (*B) with 
the charge radius of 2.5 fm is included. For the p-wave states, 
we add the Thomas-type spin-orbit interaction, with the same 
i?o and gq as of the central part. The depth of the spin-orbit 
is set to 4.40 MeV, and that of the central part is determined 
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to reproduce the proton separation energy Sp = 1.30 MeV in 
the 3/2^ state. With this potential, we have a resonance state 
at e = 0.915 MeV with the width T = 0.137 MeV in the 
1/2^ state, in good agreement with the experimental values, 
i.e.,e = 0.918 ±0.011 MeV andT = 100±20keV iH. We 
include sl/2+, pl/2-, p3/2-, d3/2+, d5/2+, £5/2", and f7/2- 
waves of the p + system in the coupled-channel calcula- 
tions. 

As for the nuclear part of the distorting potential T^^X = 
p or ®B), we adopt the microscopic folding mod el ItA. l25ll 
with the Melbourne nucleon-nucleon g matrix ll26ll . Nu- 
clear densities of *B, ^^C, ^^Al, and ^°^Pb are calculated by 
Hartree-Fock (HF) method with the Gogny-DlS force I.27L 
I28I1 . The resulting microscopic proton optical potentials are 
found to reproduce, with no adjustable parameters, the elastic 
scattering cross sections for p-'^^^Pb at 65 MeV [29i] and the 
p-^^C reaction cross sections at 200-400 MeV (HopFor ^B-A 
scattering, however, it turns out that a fine tuning of the opti- 
cal potential is necessary. This can be done with replacing the 
argument of both the real and imaginary parts of V^^' as 

rB^(l + 2;)rB, (8) 

which effectively increase the range of the potential. We 
set X to 0.04 (0.03) for the ^B-^'^C (^B-^'^Al) potential at 
285 MeV/nucleon to reproduce the experimental data of the 
reaction cross section ll3lll . As for the ^B-'^"^Pb reaction at 
65 MeV/nucleon, since there is no experimental data, we cal- 
culate the reaction cross section by CDCC with a p + ^Be + 
^•^^^Pb three-body model, and x = 0.10 is obtained to repro- 
duce the calculated value. The prescription of Eq. (|8} can be 
understood as a modification of the HF density of ^B to in- 
clude a halo structure effectively. 

TABLE II. Model space of the present calculation. See the text for 
details. 



Reaction 


elastic breakup 


proton removal 


fcmax [fm""^] 


I.O 


1.2 


Ak [fm-^] 


0.05 


0.10 


rmax [fm] 


150 


150 


-Rmax [fm] 


250 


30 


^max 


2,000 


450 



The model space of the present CDCC calculation is sum- 
marized in Table im where fcmax (?'max) is the maximum value 
of the relative wave number k (coordinate r) between p and 
®B, and Afc represents the width of the momentum bin. i?,nax 
and Lniax are, respectively, the maximum values of the rel- 
ative coordinate and the orbital angular momentum between 
^C and A. We have confirmed with the model space the 
convergence of the elastic breakup cross section (Fig. |2) for 
£ < 1 MeV and a-p (Table|IIB both within 1%. 

Results and discussion. First, we analyze the elastic 
breakup "^^^Ph^Cp^Bf^^Ph at 65 MeV/nucleon. In Fig. |2l 
we show the breakup cross section as a function of the relative 
energy e between p and ^B. We have included the experimen- 
tal efficiency e{e) ll32ll and resolution F in the calculation. We 




e[MeV] 



FIG. 2. (Color online) Breakup spectrum of the 208pb(''C,p*B)2°*Pb 
at 65 MeV/nucleon as a function of the relative energy e between 
p and *B. The dashed line shows the result of calculation with a 
normalized p-*B wave function, whereas the solid line is the result 
multiplied by I.I to fit the experimental data (J]. 



adoptr — 0.23 MeV extractedfromtheexperimentalbreakup 
spectrum of ^^C(^C,p^B)'^'^C at 65 MeV/nucleon (3^. In or- 
der to determine C^^g, we fit the theoretical result (dashed 
line) to the experimental data [4], and the solid line is ob- 
tained. The renormalization factor is 1.10, which results in 
(Cp8^B)2 = 1.78 fm-i and 5i8(0) = 67.3 eVb. 

In Fig. |2] our calculation describes well the breakup spec- 
trum below e ^ 1.0 MeV, i.e., both the transition to the 1/2^ 
resonant state and breakup to low-energy nonresonant states 
of ^C. It should be noted that we treat the resonant and non- 
resonant breakup continua on the same footing in the CDCC 
calculation. In the higher e region than the resonance energy, 
however, the calculation significantly underestimates the ex- 
perimental data. It is expected that this is due to incomplete- 
ness of our present framework. The back-coupling effects of 
three-body breakup states of^Ctop + p + ^Be on the p + ^B 
state observed will become important as e increases. In ad- 
dition, more accurate description of the p + ^B continua for 
higher partial waves with a proper p-^B interaction V^^g ^ will 
be needed. At low e, these possible problems will not exist, 
because only the tail of the overlap between ^C and p-^B con- 
tributes to the breakup process. For more detailed discussion 
on this point, see Ref. ITJ. 

To examine the peripherality of the 2°®Pb(^C,p®B)20^Pb re- 
action, we see the dependence of C^sg on the parameters of 

Vp^g ; each of Rq and is changed by 20%. Note that we put 
a constraint on the depth of the central potential so that it must 
reproduce the proton separation energy Sp. It is found that 

the uncertainty of C^fg regarding V^'g is 8%. This indicates 
that the present elastic breakup reaction proceeds peripherally 
with respect to r, as required by the ANC method. 

Second, we analyze the one-proton removal reaction of ^C 
at 285 MeV/nucleon on ^^C and ^^Al targets. As already men- 
tioned, we neglect the Coulomb breakup of ^C in this case. 
We calculate (Tbu by CDCC and the stripping cross section 
(Tstr by ERT, and obtain the one-proton removal cross section 
(T_p, as the sum of the two. Then we renormalize the calcu- 
lated (T_p to fit the experimental value taken from Ref. lIsTll . 
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TABLE ni. Results of the one-proton removal reactions with C and 
^^Al targets. The experimental data of a-p are taken from Ref. llsill . 



Target 




^■'Al 




calc. expt. 


calc. expt. 


(Jbu [mb] 
o-str [mb] 
(j_p [mb] 


2.7 
42.2 

44.9 48(8) 


4.7 
49.2 

53.9 55(11) 


{C%r [fm-i] 
5*18 (0) [eVb] 


1.73 
65.2 


1.65 
62.2 



which determines (C^Fg)^ and hence 5'i8(0). These values 
are summarized in Tal)le mil One sees that the two results of 
518(0), corresponding to ^^C and ^^Al targets, agree well with 
each other. By taking an average of the two values, we ob- 
tain (CpS^B)^ = 1.69 fm-i and 5i8(0) = 63.7 eVb. In order 
to evaluate the uncertainty of the ANC for the one-proton re- 
moval reactions, we take the same procedure as in the analysis 
of the elastic breakup reaction; the uncertainty turns out to be 
20%. By adding the aforementioned 5% uncertainty due to 
the neglect of Coulomb breakup, we find the total uncertainty 
of 518(0) extracted from a-p to be 21%. 

We here remark that in our three-body coupled-channel 
analysis, the values of 5i8(0) extracted from two different 
breakup reactions, 67.3 eVb (elastic breakup) and 63.7 eVb 
(proton removal), show very good agreement. This indicates 
reliability of the present analysis and the result of 5*18(0). As 
a principal result of the present study, we obtain (C^^g)^ — 
1.7 ± 0.3 fm^^, which corresponds to 

518(0) = 66± lOeVb. (9) 

In Fig. [3] the Sis (0) extracted by the present work is com- 

lOOi 



80- 



> 



Co 60- 



77±15 



40 



66±10 



45±13 




46±6 



• This Work 

♦ Beaumel et al. 
A Trache et al. 

X Motobayashi 

■ Descouvemont 



FTG. 3. (Color online) 5'i8(0) extracted by this work (circle) 
is compared with the results of the Coulomb dissociation method 
(cross) 0] and the analysis of a-p with the extended Glauber model 
(triangle) (3|]. Theoretical results with a cluster model calculation 
(squares) (sll and the value extracted from the d(*B, ^C)n reaction 
(diamond) [3] are also shown. 



pared with previous values. As mentioned above, previous 
results can be categorized into two, i.e., one is around 80 eVb 
(Ref. |4, 5]) and the other is around 45 eVb (Ref. i^H). Our 
result exists in between them, slightly favoring the former. 



In Ref. fT], the El contribution to the elastic breakup of 
by 208pb at 65 MeV/nucleon was extracted by subtracting 
the contributions of the nuclear and E2 breakup processes (~ 
10%) from the measured cross section, with a help of the 
breakup data by ^^C at the same energy. The rather good con- 
sistency between the present and previous results of 5i8(0) 
will indicate that the procedure for extracting the El contri- 
bution worked quite well. It was reported in Ref. [4], how- 
ever, that about 80% of the peak in the ^°^Ph(^C,p^Bf°^Ph 
breakup spectrum around e — 0.9 MeV was explained by 
nonresonant El breakup processes. On the other hand, in 
the present analysis, the peak is found to be mainly gener- 
ated by the nuclear and E2 transition to the 1/2^ resonance 
state. Reason for this large discrepancy in the resonant part 
between the present and previous studies needs further inves- 
tigation; this is our important future work. If we adopt a one- 
step calculation including nuclear and Coulomb breakup with 
all multipolarities, 5i8(0) — 54 eVb is obtained, i.e., 20% 
difference appears. This behavior is the same as in the study 
of 5i7(0) for the '^Be{p,^fB reaction (01. 

Our result is quite larger than the result of Ref. lUl, 
in which the one-proton removal reactions (^C,^B) at 
285 MeV/nucleon were analyzed by the extended Glauber 
model, with carefully evaluating the uncertainty regarding the 
nucleon-nucleon effective interactions (profile functions). By 
a detailed analysis, it is found that the difference between the 
5i8(0) obtained in the present work and Ref. |i3|] is mainly 
due to the proton optical potential. In Fig. 4 of Ref. (HI], the 
reaction cross section ctr of the p-^^C (solid line) is compared 
with experimental data. As shown in the figure, the data have 
quite large uncertainty; there seem to be two data groups be- 
tween 250 MeV and 600 MeV. Our microscopic calculation 
based on the Melbourne g matrix gives ctr = 198 mb at 
285 MeV, which is smaller than the value used in the pre- 
vious study by about 10%. It should be noted that both the 
theoretical values of ctr are consistent with the experimen- 
tal data, within their uncertainty mentioned above. This 10% 
difference is indeed crucial for the evaluation of cr_p, which 
eventually gives the difference in 5i8(0) by about 35%. Thus, 
more accurate and reliable data of ctr is highly desirable to 
judge the microscopic theoretical calculations of ctr, although 
we have shown in this study a very good agreement between 
the two 5i8(0) extracted from different breakup reactions. 

Very recently, ANCs for light nuclei with mass numbers 
between 3 and 9 are systematically evaluated by a varia- 
tional Monte Carlo calculation [33]. The resulting value of 
{CpS^}^ to be compared with ours (1.7 ± 0.3 fm^^) reads 
1.36 ± 0.03 fm^^. It will be interesting to investigate the 
difference between the two values in more detail. 

Summary. We have analyzed the elastic breakup of ^C by 
^°^Pb at 65 MeV/nucleon and the one-proton removal reaction 
of ^C at 285 MeV/nucleon on C and Al targets by a three- 
body coupled-channel framework, i.e., CDCC for the elastic 
breakup process and ERT for the stripping process. We de- 
termined the ANC C^gg and obtained the astrophysical factor 
at zero energy, 5i8(0), for the ^B(p,j)^C reaction. Our prin- 
cipal result is 5i8(0) = 66 ± 10 eVb. We have confirmed 
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that the results of S'i8(0) extracted from the two independent 
experimental data agree very well with each other, and thus 
resolved a significant discrepancy of S'i8(0) in the previous 
studies. Although the ANC is determined well in the present 
analysis, description of the breakup spectrum at higher p-^B 
relative energies is not sufficient. Extension of the present re- 
action model to incorporate the p + p + '^Be configuration will 
be very important for deeper understanding of the breakup of 
^C. Investigation on the d{^B, ^G)n transfer reaction, which 



gives a quite smaller 5*18 (0) than in the present study, will also 
be important. 
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